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ABSTRACT 

The origin of a recently discovered close-in Neptune-mass planet around GJ436 
poses a challenge to the current theories of planet formation. Based on the sequential 
accretion hypothesis and the standard theory of gap formation and orbital migration, 
we show that around M dwarf stars, close-in Neptune-mass ice-giant planets may be 
relatively common, while close-in Jupiter-mass gas-giant planets are relatively rare. 
The mass distribution of close-in planets generally has two peaks at about Neptune 
mass and Jupiter mass. The lower-mass peak takes the maximum frequency for M 
dwarfs. Around more massive solar- type stars (G dwarfs), the higher-mass peak is 
much more pronounced. These are because planets tend to undergo type II migration 
after fully accreting gas around G dwarfs while they tend to migrate faster than gas 
accretion around M stars. Close-in Neptune-mass planets may also exist around G 
dwarfs, though they tend to be mostly composed of silicates and iron cores and their 
frequency is expected to be much smaller than that of Neptune-mass planets around 
M dwarfs and that of gas giants around G dwarfs. We also show that the conditions 
for planets' migration due to their tidal interaction with the disk and the stellar-mass 
dependence in the disk-mass distribution can be calibrated by the mass distribution of 
short-period planets around host stars with various masses. 

Subject headings: planetary systems: formation - solar system: formation - stars: 
statics 
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1. Introduction 

In an attempt to place quantitative constraints on models of planet formation, we developed 
an algorithm to simulate the kinematic properties of gas giants formed in isolation (Ida & Lin 
2004a, hereafter Paper I). This prescription is based on the sequential accretion model in which 
we assume that Jupiter-mass gas-giant planets formed through 1) grain condensation, 2) runaway 
planetesimal coagulation (Greenberg et al. 1978; Wetherill & Stewart 1989; Aarseth et al. 1993; 
Kokubo & Ida 1996), 3) oligarchic growth of protoplanetary embryos (Kokubo & Ida 1998, 2000), 
and 4) gas accretion onto solid cores (embryos) (Mizuno 1980; Bodenheimer & Pollack 1986; Pollack 
et al. 1996; Ikoma et al. 2000). Based on a distribution of 1) dust-disk masses (Mj) inferred from 
mm data (Beckwith & Sargent 1996), 2) a range (1-10 Myr) of disk depletion time scale (Td cp ) 
inferred from the observed decline in the IR (Haisch et al. 2001) and mm data (Wyatt et al. 2003), 
and 3) three different growth-termination criteria (Lin &: Papaloizou 1993; Bryden et al. 1999), we 
simulated a distribution of protoplanetary masses. In addition, we considered the effect of post 
formation type-II orbital migration due to planet-disk interaction (Lin & Papaloizou 1985) which is 
an important process in relocating protoplanets away from their birth places. The last process has 
been invoked (Lin et al. 1996) to account for the origin of a population of Jupiter-mass short-period 
planets such as 51 Peg b (Mayor & Queloz 1995). 

We compared the results of our simulation with the available data of extrasolar planets. We 
suggested that around solar-type stars, 1) there may be a deficit of intermediate-mass (~ 20 — 
100M®) and intermediate-period (~ 0.1 — lyr) in their mass-period distribution (Paper I), 2) the 
frequency of gas giant planets may be an increasing function of their host stars' metallicity [Fe/H]* 
(Ida & Lin 2004b, hereafter Paper II), and 3) a large fraction of the planets migrated to the 
proximity of their host stars may have perished (Paper I, II). The first conclusion results directly 
from the expectation that a) the growth of protoplanetary cores is limited by dynamical isolation in 
the inner regions of planetary systems and slow coagulation rate in the outer regions, b) even under 
the most favorable locations, 0.1 — 1 Myr is needed for the formation of protoplanetary embryos 
(cores), anywhere in the disk, with masses M c > M ffi , c) cores with mass larger than M C;llcc ~ 
several M@ can undergo runaway gas accretion, and d) orbital migration occurs on a similar time 
scale (r m i g ) to gas depletion in the disk (rd ep ). 

The model we have presented so far provides the first step in the construction of deterministic 
properties of planetary formation. Its simplistic assumptions must be re-examined in a wider 
context of protoplanetary environment. For example, planets are assumed to form independently 
and in isolation in the model. Dynamical isolation in this context means that there are no other 
major planets around the same host stars. Most importantly, post-formation dynamical interaction 
between the planets have not yet been considered. While these effects will be examined in a future 
investigation, we consider here the mass function of dynamically-isolated planets and its dependence 
on the mass of the host star. 

On the observational side, most of the known planets are found around solar-type stars (G 
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dwarf stars) because the most successful radial velocity surveys have been conducted for these 
target stars which have cool and well defined atmospheric spectroscopic features. However, the 
search window is rapidly expanding to lower-mass stars as detection techniques are being refined 
in both radial velocity and transit searches. Figure 1 shows the distributions of semimajor axis 
(a) and mass (M p sini) of discovered extrasolar planets around M, K, G, and F stars. The planets 
around subgiants are excluded, because the relation between stellar spectral-type and its mass is 
different from that for main sequence stars. The planets discovered by transit survey are also 
excluded, because the transit survey has different observational bias in planetary periods from that 
of doppler survey. Although much more planets have been discovered around G stars than around 
F and K stars, the detection probability (after correction of metallicity dependence) is similar 
among these stars (Fischer & Valenti 2005). However, the detection probability may be by order 
of magnitude lower around M stars. Since their signals are the most conspicuous, close-in planets 
are expected to be the first to be uncovered (Narayan et al. 2004). However, until recently, the 
only planets discovered around M stars are two Jupiter mass planets with moderate semimajor 
axis (0.13AU, 0.21AU) around Gliese 876. No Jupiter-mass close-in planets have been discovered 
around M stars. 

Recently, a short-period (2.6 days) Neptune-mass (21M&) planet is found to be orbiting around 
an M-dwarf (M* = O.4M ) GJ436 (Butler et al. 2004). With a mass between those of gas giants 
and the Earth, this finding signifies a transition in the quest to search for terrestrial planets. In 
the solar system, two ice giants, Uranus and Neptune, have masses in this range. These planets 
are primarily composed of icy cores with a modest gaseous envelope. In accordance with the core 
accretion scenario, e.g., (Wuchterl et al. 2000), the cores in the outer solar system take a long time 
to emerge and when they finally acquired M C;acc , the solar nebula was already so severely depleted 
that they can only accreted a small amount of gaseous envelope (Hayashi et al. 1985). It is natural 
to extrapolate that GJ436b may also have attained M Cjacc but failed to accrete much gas. The 
main challenge to such a scenario is to account for the origin of both its low mass and short period. 

The period of the first extrasolar gas giant planet (51 Peg b) discovered (Mayor & Queloz 
1995) around a main sequence star outside the solar system is comparable to that of GJ436b. That 
planet and dozens others like it are thought to have formed through sequential accretion beyond 
the ice boundary and migrated to their present locations (Lin et al. 1996) as a consequence of their 
tidal interaction with their nascent gaseous disks (Goldreich & Tremaine 1980; Lin <fe Papaloizou 
1985). But, the much lower mass M p of GJ436b implies that gas accretion onto it may have been 
greatly suppressed prior to, during, and after its migration. 

Two other short-period Neptune-mass planets have been found around G dwarfs (a planet 
with 2.8 day period and 15M® mass around 55 Cnc and that with 9.6 day period and 14M e mass 
around HD 160691) (McArthur et al. 2004; Santos et al. 2004). Since 55Cnc and HD160691 have 
three and two other giant planets (maybe gas giants), these Neptune- mass planets could form in 
situ by accumulation of rocky materials caused by sweeping mean motion resonance associated with 
migration of a giant planet, e.g., (Malhotra 1993; Ida et al. 2000) or sweeping secular resonance 
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associated with disk gas depletion, e.g., (Ward 1981; Nagasawa et al. 2000; Lin et al. 2004). Since 
we do not include interaction from other major planets, formation of these Neptune-mass planets 
is beyond scope of the present paper. On the other hand, no gas giant planets has been found 
around GJ436. GJ436b is dynamically-isolated, so that its formation must be considered without 
help of other giant planet (s) and our calculation in the present paper can address its formation. 

Although in the present paper, we will propose a scenario that type II migration of a Neptune- 
mass planet occurs without significant gas accretion onto the planet around an M star, there 
are several other potential scenarios for the origins of dynamically-isolated, close-in Neptune-mass 
planets. Boss et al. (2002) suggested that, along with Jupiter and Saturn, Uranus and Neptune 
were formed through gravitational instability, in a massive disk. Heavy elements settled to form 
the cores and the gas envelope (clump) was greatly depleted by the photo-evaporation due to the 
UV flux from nearby OB stars before it contracts to planetary size. Since gravitational instability 
is unlikely at the present location (0.028AU) of GJ436b, the planet must have migrated from outer 
region by tidal interaction of the disk. However, the gravitational potential of disks around M 
stars is shallower than that around G stars such that the evaporation of the envelope of collapsing 
clumps would also eliminate all the residual gas in the disk. It would be difficult for the collapsing 
fragment to lose a large fraction of their mass and migrate extensively to form a GJ436b-like planet. 
On the other hand, under some extreme circumstance, the sporadic UV and X-ray irradiation from 
its host star could evaporate envelope of a gas giant planet that has migrated from outer region 
during M star's main sequence lifetime. But the overall impact of the photoevaporation process on 
planet's envelope and mass has not been determined. 

Under the general concept of sequential accretion scenario, dynamically-isolated close-in plan- 
ets may also form in situ (Bodenheimer et al. 2000). This scenario requires a concentration of 
planetesimals in the stellar proximity. One possible mechanism which may lead to such a situation 
is through embryo-disk interaction, or commonly known as type-I migration (Ward 1986, 1997a). 
The accumulation of building blocks for close-in planets also requires the termination of their mi- 
gration process and the interaction between the embryos with their host stars as well as residual 
planetesimals. 

Through a case study, we consider in this paper the origin of dynamically-isolated close-in 
planets of ~ Neptune mass around stars with various masses. With the model developed in Paper 
I and II, we show that Neptune-mass close-in planets may be abundant around M stars. In §2, we 
briefly recapitulate the sequential accretion hypothesis. Using our prescription, we first compare 
the mass function of dynamically-isolated close-in versus modest to long-period planets around 
solar-type stars (G dwarfs) in §3. This comparison is useful because we have most observational 
data and constraints for planets around solar-type stars at the present moment. Based on the best 
available observed properties of pre-main sequence evolutionary tracks of different-mass stars and 
accretion rates onto their host stars, we construct disk models around main sequence stars with 
various masses. In §4, we construct a conventional model of giant planet formation around lowest- 
mass stars and show that icy planets with 10-20M® accrete from planetesimals at ~ 1AU without 
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any significant gas accretion onto the planets. These Neptune-mass planets can also readily migrate 
to the proximity of the stellar surface, where gas accretion is quenched at Neptune-mass. Through 
a series of simulations, we show, in §5, that the disk mass dependence on their host stars' mass and 
the criteria of tidally induced migration in the disk may be quantitatively constrained by the mass 
distribution of short-period planets. We also place constraints, in §6, on the dependence of disk 
mass on the stellar mass. Finally, in §7, we summarize our results and discuss their implications. 



2. Core accretion in disks around low-mass stars 

The detailed description of sequential accretion scenario and our prescription to simulate the 
formation of planets are given in Paper I. We briefly recapitulate the central features of our approach 
and define various quantities which are used in the discussions of our results. 



2.1. Growth of protoplanetary embryos 

In the sequential accretion scenario, planetesimals grow into protoplanetary embryos (cores) 
which affects the velocity dispersion a of nearby planetesimals and modifies their own growth (Ida 
& Makino 1993; Aarseth et al. 1993; Rafikov 2003). But, a is also affected by gas drag (Kokubo & 
Ida 2002). In a disk with a surface density of dust (£d) and gas (£ g ) around a host star with mass 
M#, protoplanetary embryos' mass at any location a and time t is 

^ ( t V ( S d V ( S g \ 6/5 / m \' 2/5 / a \-9/5 fM^ 1 ' 2 

cU ~Vo.48MyrJ ) \2A x 10 3 gcm- 2 y Vl0 22 gj VlAU/ \M Q ) ® 

(1) 

where m is the typical mass of the planetesimals accreted by the embryos. Since accretion time 
Tc.grow = M c /M c increases with M c , M c (t) does not depend on its initial value M Cj o as long as 
M c (t) > M C)0 . 

In the limit of small a, the full width of embryos' feeding zone (Aa c ) is limited to ~ 10rn 
(Lissauer 1987; Kokubo & Ida 1998) where rjj is the embryos' Hill's radius rn (= (M c /3M^) 1//3 a). 
When all the residual planetesimals in an embryo's feeding zone have coagulated with it, the embryo 
attains an isolation mass 

— -Kia-) 3/2 (^) 3 (^) 3/2 (t)" 1% »- 



2.2. Protostellar disk properties 



Equations (1) and (2) indicate that both M c (t) and M C) i so are determined by the distribution 
of Ed and S g . In Paper I and II, we introduced a multiplicative factor (/d and / g ) to globally scale 
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disks with the minimum mass nebula model for the solar system (Hayashi 1981) such that 



where the step function r][ cc = 1 inside the ice boundary at aj ce and 4.2 for a > aj cc . [Note that the 
latter can be slightly smaller (~ 3.0) (Pollack et al. 1994).] The minimum mass model corresponds 
to /d = / g ~ 1. Here we introduced a new scaling factor, /id and h g , representing the dependence 
on stellar mass (see below). 

If the disk is optically thin and heated by stellar irradiation only (Hayashi 1981), 



The stars' luminosity L* is generally a function of their mass M* and age t. Additional heating due 
to viscous dissipation enlarges this boundary (Lin & Papaloizou 1980). This ad hoc phenomeno- 
logical prescription provides a useful working hypothesis for comparative analysis between solar 
system architecture and extrasolar planetary systems. Since is generally an increasing function 
of M*, CLi ce is small in disks around low-mass stars. 

Around solar-type T Tauri stars, the observationally inferred total mass of dust, Md, in the 
protostellar disks ranges from 1O _5 M to 10 _3 Mq (Beckwith &: Sargent 1996), which corresponds 
to a range of /d ~ 0.1 — 10. The disk-mass determinations are somewhat uncertain due to the 
poorly known radiative properties of the grains. Nevertheless, their divergent dust content provides 
a reasonable evidence for a greater than an orders of magnitude dispersion in M^. In addition, 
the radio image of the disks is not well resolved in many cases. A rough magnitude of /d can be 
inferred from the the total mass of the disk under the assumption that all disks have similar sizes 
(a few tens AU to a hundred AU). In this paper we follow Paper II and generate a set of /d with a 
unit variant of Gaussian logarithmic distributions and a range between 0.1-10 with cut-off of high 
fdhd tails at 3O(M*/M ) (see Figure 4). If f g h g > 30(M*/M & ), the disk is gravitationally unstable 
at a > a few AU (see discussion in section 7), and the total mass of heavy elements in the disk 
would be significant fraction of that in the host star. 

In eq. (3), the dependence of Ed on M* is incorporated in a mass scaling function h<±. In 
general, a relatively weak IR excess is associated with low-mass stars, which suggests /id may be an 
increasing function of M*, at least in the inner regions of the disk. For host stars with M* = M Q , 
h d = l. 

Equation (1) indicates that M c is also a function of S g through damping of a due to gas 
drag. The planetary migration rate that we mention later is also dependent on S g . Around other 
stars, very little information is available on S g . Since there is no indication of divergent depletion 
pattern between molecular hydrogen and mm-size dust emission (Thi et al. 2001), we follow the 
same prescription for Sd with a disk mass scaling parameter f g and a stellar-mass dependence 
function h g (eq. [3]). Similar to Paper I and II, we adopt the conjecture that fd does not change 





= 2.7(L*/L ) 



V 2 AU. 



(4) 



- 7- 



except in those regions where they have been totally accreted by the cores and 

ig = / g ,oexp( — J , (5) 

V r dcp/ 

where Td ep is disk depletion time scale discussed below. The assumption of the uniform exponential 
decay is for simplicity. The effects of detailed decay pattern as a result of viscous evolution of disks 
will be discussed in a separate paper. For computational simplicity, we assume [Fe/H]=0, i.e., a 
solar composition for all stars so that /d = / gi o and h g = h^. The dependence of planetary M p — a 
distribution on [Fe/H] for solar- type stars has already been discussed in Paper II. 

We now consider a prescription for the stellar mass dependence. The parameters are h& 
(= h g ) and r^ cp . In young clusters, the fraction of stars with detectable IR (Haisch et al. 2001) 
and mm continuum ( Wyatt et al. 2003) from circumstellar disks around T Tauri stars declines on 
the time scales of 1-10 Myr. Although this decline may be due to dust growth and planetesimal 
formation rather than the depletion of heavy elements (D'Alessio et al. 2001; Tanaka et al. 2005), 
the correlation between the intensity of mm dust continuum with the gas decline emission (Thi 
et al. 2001) and the UV veiling for ongoing gas accretion suggest that gas is depleted as the dust 
signature fades. The dust signature maintains up to lOMyr also for disks around brown dwarfs 
(Mohanty et al. 2003). So, in this paper, we assume Td ep in the range of 1-10 Myr for all stellar 
masses. 

Although direct estimates of disk mass both in gas and dust are difficult to obtain for the 
inadequate sensitivity of existing observational instruments, the disk accretion rate M can be 
inferred from the Ha line profiles (Muzerolle et al. 2003; Natta et al. 2004) such that M oc M 2 with 
a large dispersion. If the angular momentum transfer and mass diffusion time scale is insensitive 
to M*, we could infer h g (= /id) ~ (M*/M Q ) 2 . In view of the large uncertainty in the data, we 
consider three possible dependence on M*, 

h d = (M./Mq) ' 1 ' 2 , (6) 

with hd = (M^/Mq) 2 as a standard case. 



2.3. Core growth, isolation mass, and gas giant formation around stars with 

different masses 

With these prescriptions for disk parameters, we find that 

by assuming m = 10 22 g. From eq. (3), we estimate the core masses in protoplanetary systems with 
cq. (2) such that 

Me, o ,0, 6 ^/f*f(^) 3/1 (t)" 1/2 Me. <S) 
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We have already pointed out in Paper I and II that the growth of planetesimals is limited by 
isolation at small a and the slow coagulation rate at large a. Since for smaller M#, rn is larger, 
M Ci ; so is larger for the same mass disks. However, both hd and h g are increasing functions of M*, 
so that the growth of M c (t) and the isolation mass M Cj j so at any given a are actually smaller for 
lower-mass stars. We do not explicitly include type I migration (Ward 1986, 1997a) of cores, but 
take into account its effects in some runs (see Paper I and discussion in §5). 



As the cores grow beyond a mass 

0.25 

'c, hydro 



M Cihydro ~ 10 [ — ^ — ] M e , (9) 



W- 6 M (B /yr / 

their planetary atmosphere is no longer in hydrodynamic equilibrium and they begin to accrete 
gas (Stevenson 1982; Ikoma et al. 2000). In the above equation, we neglected the dependence on 
opacity (see Paper I). In regions where they have already attained isolation, the cores' accretion 
M c is much diminished and M Ci h y dro can be comparable to an Earth mass. But the gas accretion 
rate is still regulated by the efficiency of radiative transfer such that 

IMzjl^Me. no) 

dt TKH 

where M v is the planet mass including gas envelope and the Kelvin Helmholtz contraction time 
scale is (for details, see Paper II) 

™ ^ 10 '° OS)" ^ (11> 

Gas accretion onto the core is quenched when the disk is depleted either locally or globally. 
A protoplanet induces the opening of a gap when its rate of tidally induced angular momentum 
exchange with the disk exceeds that of the disk's intrinsic viscous transport (Lin h Papaloizou 
1985), that is, when the planet mass M p exceeds 

M 5 ,™ , £ M, ^ (i ) 2 M„ 3 (^) (^) * (£) V1 M e , (12 ) 

where we used an equilibrium temperature in optically thin disks (Hayashi 1981) and a-prescription 
for the effective viscosity v (Shakura & Sunyaev 1973), in which T = 170(a/ai cc ) -1 / 2 K and 
v = ah 2 ^K where a is a dimensionless parameter, h and f^K are the disk scale height and Kepler 
frequency. Since h = c s /Qk where c s is the sound velocity, h 2 oc T/M* oc L l J 4 /M*, so that 
Mg iV i s oc L l J A . Since is an increasing function of M*, eq. (12) indicates that M gjV i s is smaller 
around lower-mass stars. The orbital evolution of planets is locked to the viscous evolution of the 
disk gas (type II migration) when their M p > M Pim i g = A v M g ^ s where A v is the dimensionless 
factor ~ 3 — 10 for a laminar disk (Lin & Papaloizou 1985). 
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The type II migration rate a is given by a/r m i g with 

Tm ^ 10 v E -v(^)(t)(ilu) 1/2 - (») 

where Mj is Jupiter mass. We set a lower limit on r m i g at a 2 /V ~ 4.3 x 10 5 (a/10~ 4 )(a/lAU) yrs 
(Paper I). We found that the formula for J m (angular momentum flux at the radius of maximum 
viscous couple, r m ) given in Paper I (eq. [63] in the paper) must be multiplied by a factor 2tt. If 
we use J m to evaluate evolution of planetary orbital radius a, eq. (13) is reduced by 2ir as well 
(Paper I). However, since the planetary migration may be caused by a fraction of J m and the 
fraction is uncertain, we use eq. (13) in the present paper too. We adjust the migration time 
scale by the value of a, which is also uncertain, comparing with observational data. As shown in 
Paper II, in order to reproduce period distribution of gas giant planets similar to that of observed 
extrasolar planets around solar-type stars, the disk viscous diffusion time scale r v = r^/v ~ 
4 x 10 6 (a/10~ 4 )(r m /10AU) yrs must be comparable to disk lifetime Td ep . Hence, the a viscosity 
in our model must be ~ 1CT 4 . Although the regions at 1-10AU could be "dead zone" for MHD 
turbulence (Sano et al. 2000), resulting in a very small a in these regions, the best fit value a ~ 10~ 4 
may not necessarily reflect a realistic value because of rather simple assumptions for E ff distribution 
and its exponential decay in our model. We will carry out more detailed calculations coupled with 
disk viscous evolution in a separate paper. 

The planets' migration is terminated either when the disk is severely depleted (/ g — ► 0) or when 
they reach a sta ii, which is set to be 0.04AU in our calculations. There are potential mechanisms to 
stop migration at < 0.05AU (Lin et al. 1996). However, Paper II suggests that only a small faction 
(< 10%) of migrating planets can survive in the vicinity of their host stars. This fact should be 
kept in mind when population of close-in planets is discussed with our model. 

The gap becomes locally severely depleted when the planets' Hill's radius (rn) exceeds the disk 
thickness h (Bryden et al. 1999), that is, when M p exceeds 

Growth through gas accretion is quenched for planets with M p > M P) t runc = ^th-Mg,th- Numerical 
simulations show some uncertainties in the dimensionless parameter A t ^ (Bryden et al. 1999; Nelson 
et al. 2000). Planets with M Pitrunc > M p > M PjI11 i g migrate with the disk while continue to accrete 
gas, albeit at a reduced rate, e.g., (Lubow et al. 1999). We use eq. (10) without a reduction factor 
for simplicity, because the reduction factor is quite uncertain and introduction of the factor does 
not affect the results significantly. Since gas accretion for Mp > -^p,mig 

is already very rapid (11), 

the reduction does not change total gas accretion time scale. Following Paper I, we adopt in this 
paper A t h = 1.5 3 ~ 3.4, that is, the truncation condition is rn > 1.5h. Since rapidly increases 
with M*, M P)trunc is also smaller around lower-mass stars. 

Even for planets with M p < M Pi t rU nc) gas accretion may be ultimately limited by the dimin- 
ishing amount of residual gas in the entire disk. For our disk models, the maximum available mass 
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M g)noiso ~ vra 2 S g ~ 290f gfi h g ^ exp M ®- (15) 

When M g noiso becomes smaller than M p , gas accretion is terminated. 

A similar global limit M C)I10 i S0 = 7ra 2 £d (see eq. [20]) is also imposed if the formula (2) exceeds 
it. As the gas is severely depleted, the velocity dispersion a of the embryos and residual plan- 
etesimals grows until they cross each other's orbits (Iwasaki et al. 2002; Kominami & Ida 2002). 
Eventually a few surviving embryos acquire most of the residual planetesimals and less massive 
cores during the late oligarchic-growth stage. The asymptotic embryos' masses are given by eq. (2) 
with Aa c ~ Kmrf/^K, 

where V SUT f and pd are the surface escape speed and internal density of the embryo. We use this 
enlarge asymptotic mass when f g < 10~ 3 . 

In Paper I and II, we put all of these processes into a numerical scheme to simulate the 
formation and migration probabilities of planets around solar- type stars. Cores with M p > M C]3CC 
emerge on time scale shorter than Td cp in disks with modest-to-large values of /d around solar-type 
stars. In this limit, gas accretion and orbital migration lead to the formation of gas giants with 
kinematic properties similar to those observed. 



3. Emergence and migration of Neptune-mass planets in disks around solar-type 

stars 

We now apply our numerical methods to study the formation of planets around stars with 
various masses. Our objective is to simulate the mass function of close-in planets and assess the 
influence of formation and migration on it. We show here that this quantity can provide clues on the 
dominant processes which regulate planet formation and it also can be used to distinguish between 
competing theories of planet formation. In all models, we choose a = 10~ 4 based on assumption 
that the viscous evolution time scale for the disks is comparable to Td ep (see §2.3). Using these 
models, we carry out Monte Carlo simulations. For simplicity, we generate a set of initial a's of 
the protoplanets and r^ cp with uniform distributions in log scale in the ranges of 0.1-100AU and 
10 6 -10 7 yrs. The assumed a distribution corresponds to orbital separations Aa that is proportional 
to a, which may be the simplest choice. The distribution of /d was discussed in §2.2 (also see Figure 
4). In Paper I and II, we also assumed a distribution of M* in a range of O.7-1.4M . In this paper, 
in order to make the M* dependence clear, the value of M* is fixed in each run. 

In the following sections, we consider several sets of model parameters. We first discuss a 
standard model with A v = 10 and = (M*/M ) 2 around a solar-type star with M* = l.OM 
(model 1.0). The evolution of totally 20,000 planets are calculated for each run. 
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3.1. Formation of cores and asymptotic mass of protoplanets. 

Model 1.0 is similar to the results we have already presented in Paper I. In Figure 2a, we 
highlight the mass (M Pinn ) and semi major axis (a nn ) distribution of planets at t = 10 9 yrs after 
they have attained their asymptotic mass and gone through the initial migration due to their tidal 
interaction with their nascent disks. The main features to notice in this panel are: 1) a deficit of 
planets with intermediate masses (20 — 100M®) at intermediate semi major axis (0.1 — 1AU), and 2) 
a large population of close-in (a nn < 0.05AU) gas giants with 2OM0 < M p <2x 10 3 M®, although 
only a small faction of them (< 10%) may be able to survive (Paper II). We have already indicated 
in Paper I that these properties are due to 1) the runaway nature of dynamical gas accretion and 
2) type II migration. 

In order to distinguish between these two dominant effects, we trace back, in Figure 2b, 
the initial semi major axis (a- m i) where the cores of both close-in (marked by black circles) and 
intermediate or long-period planets (marked by gray dots) formed. These results clearly indicate 
an one-to-one mapping between the mass function of the close-in planets and the locations where 
they are formed. For illustrative purposes, we also mark the domain where some physical processes 
operate and dominate the evolution of protoplanets. For example, the thin solid lines indicate 
the upper limit of the isolation mass that cores can attain prior to gas depletion (it is obtained 
with /d = 30 in M Cj i so and M c ^ no - mso ). The transition at 2.7 AU in model 1.0 corresponds to the 
ice boundary. The thick solid lines indicate the critical mass for the onset of type II migration, 
-^p,mig = A u M gtV i s with A v = 10 according to eq. (12). We also highlighted the asymptotic growth 
limit M P)trunc = AfoMg^h with A t ^ = 3.4 (eq. 14) by broken lines. 

At any given ai n i, a fraction of terrestrial planets can form with M p above M Pi i so through 1) 
gas accretion and 2) merger of residual planetesimals and other embryos after disk gas depletion. 
When the embryos reach their isolation mass, M c h y dro declines with the vanishing M c (eq. [9]). 
Although gas accretion is initiated, planets cannot grow significantly prior to severe gas depletion 
unless the planets' isolation mass M C) i so > M Cjacc ~ several M®. The gas accretion time scale 
given by eq. (11) can be comparable to or shorter than Td ep only for cores with M c > several M®. 
Mergers of residual planetesimals and other embryos occur during and after the gas depletion (see 
eq. [16] for asymptotic M c ~ M ej i so ) (Kominami & Ida 2002). 

If their M p > M Pjm i g = A u M g y m and there is adequate residual gas in the disk (f g is still 
large enough), they would migrate to the vicinity of their host stars (eq. [13] with eq. [5]). The 
magnitude of M giV j s is an increasing function of a. For a < 0.7 AU, tkh for embryos which can 
migrate is longer than their r m i g . Although their M p may be < M Pitr unc = A;h-^g,th> their inward 
migration is sufficiently rapid that they do not acquire a significant amount of gas along the way 
(Ivanov et al. 1999). Termination of migration allows low-mass planets (formed at a < 0.7AU in 
relatively massive disks) to grow through gas accretion to their asymptotic mass M Pjtr unc — 30M® 
at the stalling location which is set to be a s t a ii = 0.04 AU. Thus, the final mass of close-in planets 
with a; n i ^5 0.7AU is ~ 30M® regardless of truncation mass M Pitr unc at the original locations ai n i. 
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These planets are formed interior to the ice boundary and they are likely to be mostly composed 
of silicates and iron, in contrast to the ice giants in the solar system. Also, they cannot accrete 
large amount of gas because the aspect ratio of their nascent disk at 0.04 AU is so small that they 
induce clear gap formation when their mass reaches that of Neptune. 

Planets formed at slightly larger o's (~ 1AU) must attain M p > 20 — 30M® before they acquire 
the mass to start migration, M Pim i g = 10M giV i s (eq. [12]). With this critical mass, embryos which 
can initiate migration there can also accrete gas efficiently with tkh < fmig provided the gas in the 
disk is not severely depleted. Since M Pjtrunc decreases with decrease in a, when they arrive close to 
their host stars, their M p may be > M Pttrunc so that they would not acquire any additional mass. 
Note that cores of these planets are also mostly made of silicates and iron rather than ice and they 
should not be referred to as hot Neptunes. 

At even larger radius, M Ci i so > 5 — 10M®. Prior to reaching isolation, embryos' gas accretion 
is suppressed by the bombardment of residual planetesimals, i.e., M c < M c hydro f° r moderate M c . 
After isolation is reached, M Ci h y dro becomes smaller than M c and their tkh due to gas accretion is 
reduced below r m i g and they quickly evolve into gas giants. Thus, final M p in Figure 2b coincides 
with M Pi trunc at initial locations. Beyond ~ 10 AU, the time scale for the emergence of cores with 
M p > M Cjacc or isolation is comparable to or longer than both T m - lg or Td cp . Although they may 
acquire M p in the range of 10 — 100M® through merges of residual embryos after the gas depletion, 
these planets generally do not migrate extensively. 

3.2. Mass distribution of short-period planets 

Interior to the ice boundary, nearly all the planets with sufficient mass to initiate efficient gas 
accretion have migrated to the vicinity of the host star (see Figure 2b). But, a majority of the 
planets which migrated to the vicinity of solar-type stars were formed beyond the ice boundary, as 
gas giants, prior to their migration. There is a narrow window in the range of a where the seed 
of intermediate-mass planets may form and migrate to the proximity of their host stars. In Figure 
3, theoretically predicted mass distributions are plotted. Since it is expected that most of close-in 
planets may fall onto their host stars, we plot the distributions of close-in planets, reducing the 
amplitude N by a factor 10. For comparison of the amplitude between close-in and distant planets, 
uncertainty in this calibration is noted. We also plotted observed distributions. Since the number 
of runs in each model does not reflect the number of targets for current doppler survey, we cannot 
compare the amplitude N between the observed and the theoretically predicted distributions. Only 
the shape of the distributions should be compared. Also note that observed distributions do not 
exactly correspond to host stars' mass of each model and numbers of observed planets are not large 
enough for statistical arguments for stars other than G stars (model 1.0). We show that, around 
stars with M* ~ 1M Q , the mass distribution for the close-in (with ag n < 0.05 AU) planets is skewed 
toward ~ 10 3 M® (model 1.0 in Figure 3a). This distribution is more enhanced near ~ 10 3 M® than 
that observed. The effect of post formation star-planet tidal interaction, which has not been taken 
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into account in our model, may have caused the demise of a majority of the close-in planets, in 
particular massive planets (Gu el al. 2003, Paper II). 

In disks with modest masses, planets form in the advance stages of evolution when gas depletion 
is well underway Some of these planets may migrate interior to the ice boundary and become 
stalled while others remain close to their place of birth. A majority of gas giant planets with 
0.1 AU < ag n < 1 AU have migrated but not extensively. In contrast to the close-in planets, 
model 1.0 in Figure 3b clearly shows a paucity of longer-period (0.1AU < ag n < 1 AU) planets at 
M p ~ 20 — WOMq. This distribution reflects the stringent prerequisite that gas accretion into gas 
giants must be preceded by the rapid formation of sufficient mass cores whereas the build up of 
terrestrial planets can continue well after the severe depletion of the disk. 



In a generalization of the solar nebula model, Nakano (1988) showed that the temperature 
distribution throughout the disk increases with M*. However, he did not consider the dependence 
of Sd and S g on M*. In a recent paper, Laughlin, Bodenheimer & Adams (2004), considered a 
model in which the disk mass increases with the stellar mass. Their objective is to demonstrate 
the difficulties to form gas giants around M dwarf stars. But, the effects of planetary migration, 
truncation of gas accretion due to gap opening and the gradual depletion of the disk gas are 
neglected. 

In this section, we consider the variation of 3 model parameters: 1) the stellar mass M*, 2) 
the dependence of disk mass on the stellar masses, hd, and 3) the condition for the onset of type 
II migration. A standard series are Model x where x = 0.2, 0.4, 0.6, 1.0, 1.5 represents models with 
M* = xMq. In the standard series, we set A v = 10 and h& = (M*/Mq) 2 . Model xB (series B) 
and xC (series C) correspond to A v = 1 and 100 with h& = (M^/Mq) 2 . x in xB and xC expresses 
M*/M Q as well. Model xD (series D) and xE (series E) correspond to h& = M*/M Q and 1 with 
A v = 10. 

In Figure 4, distributions of f^hd we used for models x and xB are shown. They represent 
the relative mass distribution of disks around stars with various masses. The mean value of Sd 
is an increasing function of M* (see eq. [3]). In order to limit additional model parameters, we 
assume ZAMS mass-luminosity relationship, L*/Lq = (M*/M@) . Since the time scale of pre-main 
sequence stage of lower-mass stars is long, planet formation around these stars may proceed during 
their pre-main sequence stage in which L* is rather large. In that case, the dependence of 
on M* is weaker, but it may still have a positive power-law dependence, so that the trend of the 
M* -dependence of planetary systems shown below does not change. 

In the calculations in this section, 
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4.1. Planets around M dwarf stars 

In Figures 2 and 3, we also included the results of the simulations for models 0.2-0.6 and 
1.5. We first present a low-mass model 0.2 since it is in strong contrast to model 1.0. Stars with 
M* = 0.2Mq corresponds to relatively light M stars. These stars are not only most numerous but 
they also contribute most to the initial stellar mass function. According to our prescription and 
model parameters, £ g and around the host star in model 0.2 are 25 times smaller than those in 
model 1.0. 

In the result of model 0.2 on the top panels of Figures 2a and b, we find that Jupiter-mass 
planets rarely formed around low-mass stars. This paucity is due to the slow growth rate of embryos 
and their low isolation mass such that little gas can be attained by them prior to its depletion. 
This result confirms the conclusion reached earlier by Laughlin, Bodenheimer <fe Adams (2004). 

The upper limit of the M Pi fi n — a^ n distribution is determined by the M Ci ; so . This correlation 
arises because is a rapidly rising function of M*. In our prescription, the ice boundary is 
located at ai ce ~ 0.11AU (eq. [17]) which is 25 times closer to a host star with M* = O.2M than 
in model 1.0 with solar-type stars. Nearly all the cores formed around these low-mass stars are 
mostly composed of ice. Planets can acquire masses M p > M Cj i so through gas accretion. But with 
a relatively small M Ci ; so , the accretion is inefficient (eq. [11]). They can also gain mass after the 
gas depletion through collisions and mergers of residual embryos. But the asymptotic mass M ei i so 
is also relatively small for low-mass stars (eq. [21]). 

With the ad hoc a prescription we have adopted, the necessary condition for the onset of type 
II migration is satisfied for relatively low-mass planets (eq. [22]). When their M p > M P)m i g , these 
cores undergo orbital decay. Similar to model 1.0, there is a population of close-in planets with 
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2M e < M p < 10M®. The results on the top panel of Figure 2b show that they indeed originated 
from a region between ai ce and ~ 3 AU. They also indicate that planets with mass down to 2M§ 
may migrate to the proximity of a 0.2M Q star. At the arbitrary a s t a ii = 0.04 AU around such a 
host star, the equilibrium temperature of these close-in Neptune-mass planets is similar to that of 
the Earth. The composition and structure of such planets have been already been discussed by 
Leger et al. (2004). 

Although a few planets formed with M p > 10M e , they are the exceptional cases resulting 
from the tails of the fd distribution. Consequently, the mass function of planets with modest or 
large a^ n (> 0.1 AU) show a sharp decline at M p ~ 10M ffi (Figure 3b). It also extends well into 
the low-mass range. In contrast, the mass function for the close-in planets (with afi n < 0.05 AU) 
shows a peak near M p ~ 10M® (Figure 3a). 

In model 0.4, we consider a host star with M* = 0.4M Q which corresponds to a relatively 
massive M dwarf star. It is an analog of GJ436, around which a single close-in Neptune-mass 
planet has been discovered (Butler et al. 2004). Figure 3a and b indicate that the frequency of 
planets observable with current doppler survey is significantly smaller for M stars than for K, G, 
and F stars (also see Figure 8). This is consistent with the observation and the finding of Laughlin, 
Bodenheimer & Adams (2004). 

Figure 5 indicates the mean mass and characteristic mass associated with the highest peak of 
the mass distribution for a) the close-in planets and b) planets with 0.1 < a < 1AU. These masses 
increase with M*, in particular for the close-in planets. This accounts for the recent planetary 
finding around GJ436. 

We find the frequency of Neptune-mass close-in planet peaks at M* ~ 0.4M© (Figures 3a and 
5). The formation of these planets depends sensitively on the environment and their frequency 
provides constraints on the sequential accretion scenario. Most of these planets formed slightly 
beyond the ice boundary. For this stellar mass, the ice boundary is located at ~ 0.43 AU. Outside 
it, the upper limit of M Cj i so exceeds M Pjm i g . In Figure 6, we illustrate the evolution of a typical 
embryo that forms near the ice boundary in a disk with h^fd = 7.2 and a depletion time scale 
Tdep = 9.2 Myr. The disk surface density h^fd = 7.2 is close to the tail of the distribution for 
OAMq stars (Figure 4). During the initial 2 x 10 4 yrs, an embryo grows to ~ 10M© through 
coagulation. The rapid coagulation enhances M Cj h y dro) which prevents gas accretion (see eq. [9] 
and Figure 6). Since this embryo is located beyond the ice boundary and hdfd = 7.2, M Cj i so is 
as large as ~ 30M e . But, before it acquires M Cj i so , it opens up a gap and undergoes type II 
migration at M c ~ 10M e . During the migration, gravitational perturbations from the embryo 
prevents additional planetesimals from reaching the embryo and its growth is quenched (Tanaka 
& Ida 1997; Rafikov 2003). On a time scale of 3 x 10 5 yr, the embryo migrates to 0.04 AU where 
it is stalled. The termination of bombardment by residual planetesimals makes M c h y dro lower to 
start gas accretion onto the embryo. However, the gas accretion time scale tkh is longer than r m i g 
(Figure 4). Gas accretion onto the planet actually proceeds after it reaches 0.04AU. The planet's 
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growth is quenched at ~ 14M® as a clean gap is formed because M Pjtrunc ~ 14M® at 0.04AU. 
Figures 2 and 3 show that many planets evolve in a similar way to acquire Mpg n ~ 14M®. The 
composition of these planets is similar to that of Uranus and Neptune and its day-side surface 
temperature may be ~ 500K. 

In the lowest-mass model 0.2, the upper limit of M Cj i so does not exceed M Pjm i g even outside the 
ice boundary. The migration of Neptune-mass planets requires a delicate balance between M Ci i SO 
and M Pjm ig. Thus, the frequency of close-in Neptune-mass planets in model 0.4 is larger than that 
in model 0.2. 

This frequency in model 0.4 is also larger than that in model 1.0. Although the upper limit of 
M Ci ; so is larger than M Pim i g outside the ice boundary around solar-type stars, M p > 30 — 40M® at 
the onset of migration, so that gas accretion is much more efficient onto such large cores (eq. [11]). 
Thus, most isolated planets tend to arrive at the proximity of a solar-type star as gas giants and 
only a small fraction arrives as silicate and iron cores with limited gaseous envelopes. 

Recently Neptune-mass planets have been detected around 55Cnc and HD160691 (McArthur 
et al. 2004; Santos et al. 2004). Since these stars are G-type stars, if the result of Figure 3 is applied, 
the probability of formation of such planets is very low although it is not zero. Unlike GJ436, there 
are three additional Jupiter-mass planets around 55Cnc and two additional ones around HD160691. 
In the system of 55Cnc, sweeping mean motion resonance associated with migration of the giant 
planet (Malhotra 1993; Ida et al. 2000) presently at 0.1AU could bring rocky embryos/planetesimals 
to the vicinity of the host star, so that a Neptune-mass rocky planet could accrete in situ. On the 
other hand, in the system HD160691, sweeping secular resonance associated with disk gas depletion 
(Ward 1981; Nagasawa et al. 2000), could bring rocky embryos/planetesimals to inner regions (Lin 
et al. 2004). In the present paper, we do not include such interactions. We will present elsewhere 
details. In the system of GJ436, however, no additional Jupiter-mass planet has been found. Our 
model accounts for formation of the isolated Neptune-mass planet around an M star. 



4.2. Planets around K dwarfs 

Figure 3b shows that in model 0.6 where M* = O.6M , the intermediate-mass (20 — 100M®) 
and intermediate-a (0.1 — 1AU) planets are more abundant than in either model 1.0 or model 0.4. 
The mass function of close-in planets also appears to be smoother with some M p in the range of 
~ 10 — 100M® (Figure 3a). These intermediate-mass planets formed just outside the ice boundary 
in disks with modest fa where the upper limit of M c [ so ~ 50 — 100-/W® while that for typical disks 
(with fa ~ 1) is ~ M®. The critical mass for starting planetary migration at the ice boundary is 
20M® in model 0.6. A small fraction of emerged cores may accrete modest amount of gas as it 
starts to migrate. 

The formation of the intermediate- mass and intermediate-a planets, which tends to smooth 
the mass distribution, may be one of characteristics of planets around K stars, compared with 
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those around G and F stars, although it is much less pronounced than the characteristics of planets 
around M stars. For illustration, we present, in Figure 7, the formation of a typical planet which 
formed with an intermediate mass and attain an intermediate a during its migration. In this case, 
a seed embryo is formed at aj n i = 1.8AU in a slightly massive h^fd = 2 disk with Td ep = 2.3Myr. 
Through planetesimal coagulation, this embryo attains a mass M p ~ 8M e and become isolated in 
1.5 Myr. The cessation of the planetesimal bombardment enables the embryo to grow through gas 
accretion. When its mass reaches M p = 30M® at ~ 8 Myr, gas accretion is quenched by the severe 
depletion of gas near its orbit. The newly formed planet undergoes migration while gas is globally 
depleted. The orbital migration is eventually halted at an intermediate location, 0.4AU. 

The condition for growth to be quenched between 10 — IOOMq requires a ~ 1 AU. The above 
example shows that to halt migration at an intermediate location, both time scales of migration 
and growth due to gas accretion are required to be comparable to the gas depletion time scale. In 
general, such special circumstances are satisfied with small probability. However, they are more 
likely around K dwarfs than other type stars because cores with M p > 10M® are more abundant 
at a ~ 1AU around K dwarfs as a result of the M* dependence of ai ce and M Ci i co . In model 0.6, 
gas accretion may be quenched when the cores attain a mass lOM® < M p < 100M® while they 
migrated to 0.04AU < a < 1 AU. In comparison with model 1.0, the isolation mass M C) i so at the ice 
boundary increases with M*. Around a solar- type star, the isolation mass near the ice boundary is 
sufficiently large for efficient gas accretion to be initiated. Eventually runaway gas accretion leads to 
the emergence of the intermediate-mass deficit in the mass distribution of planets around relatively 
high-mass stars. Figure 1 shows that some fraction of planets discovered around K stars may 
have the intermediate mass. However, the number of the detected planets may be insufficient for 
statistical discussion. These stellar mass dependence in the extrapolated planetary characteristics 
can be tested with future observation. 



4.3. Planets around higher-mass stars 

In model 1.5 where M* = 1.5M Q , the range of a where gas giants are formed is more extended 
(the bottom panel in Figures 2). This arises primarily because the disks around more massive stars 
have relatively large Ed- In this case, M Cjacc (the core mass required for rapid gas accretion) can 
be attained before gas depletion even at large a. However, larger aj ce in this case leads to slow core 
growth beyond the ice boundary. Thus, most gas giants have cores composed of silicates and iron 
but not icy cores, in contrast with gas giants around lower mass stars. The very laxg6 dice 

also leads 

to less efficiency of formation of gas giants in the range of a few AU to 10AU than that around G 
stars, although gas giants form in broader range of a. As a result, the fraction of F stars harboring 
giant planets with periods smaller than several years which are currently detectable with doppler 
survey is rather smaller than that of G stars (Figure 8), although the difference is within a factor 
1.5. The similar fraction within a factor 2 among K, G, and F stars shown in Figure 8 is consistent 
with the observation (Fischer & Valenti 2005). The predicted mass distributions for M, K, G, and 
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F stars are not inconsistent with the observed distributions. Since the numbers of observed planets 
are not enough for statistical discussion, in particular for stars other than G stars (Figure 3), we 
cannot discuss the agreement between the predicted and observed distributions in more detail. 

Strong winds and jets from further higher mass stars (A, B stars) may decrease Td ep , which 
reduces formation rate of gas giants. It is not clear how much the fraction of massive stars with 
gas giant planets is reduced. We will address planetary formation around massive stars elsewhere. 

5. Dependence on the migration condition 

The results of the standard models clearly indicate that the mass function of planets depends 
on the delicate balance between growth and migration time scales. There are some uncertainties 
concerning the migration process. In the standard series, we set A v = 10 (eq. [22) in accordance 
with the results of previous numerical simulations (Lin & Papaloizou 1985). However, the additional 
contribution from a torque imbalance between the Lindblad resonances of low-mass embedded cores 
may also lead to type I migration while their masses are relatively small (Goldreich & Tremaine 
1980; Ward 1986, 1997a), although turbulence in the disk (Nelson & Papaloizou 2004; Laughlin, 
Steinacker & Adams 2004) and self induced secondary instability (Balmforth & Korycansky 2001; 
Koller et al. 2003) can also retard the rate of migration. In our simulations, we can partly take 
into account the effect of type I migration by lowering the value of A v (to unity). We also carried 
out simulations with A v (= 100) such that the onset of migration is delayed until the cores have 
attained relatively large masses. 

The sensitive dependence of the mass function of close-in planets on the migration condition 
makes it an ideal observable feature which can be used to calibrate the criteria and efficiency of 
migration. In a variation of the standard models, we consider two new series: 1) series B (model 
xB) with A v = 1 and 2) series C (model xC) with A v = 100, which have identical M* as the 
standard series with A v = 10. In Figures 9a and b, we show the final mass and semimajor axis 
distribution of planets in models xB and xC. 

For models xB, cores undergo migration before they attain sufficient mass to engage in efficient 
gas accretion. Although the cores' migration may terminate close to their host stars, the relatively 
small aspect ratio of the disk for small a implies low M Pj t runc . At small a, gap formation prevents 
the cores from accreting gas. The accumulation of cores and planetesimals in the proximity of 
their host star may promote their coagulation (Ward 1997b). Although we cannot rule out the 
possibility of a highly-efficient migration on the basis of the mass function of the close-in planets, it 
does pose difficulties to account for the modest frequency of gas giants with periods longer than a 
few weeks. Even around solar-type stars, cores rapidly migrate to the stellar proximity before they 
have acquired sufficient mass to efficiently accrete gas so that the probability of gas-giant formation 
is strongly suppressed. 

For the low-stellar-mass models 0.2B and 0.4B, the isolation mass is only a few times larger than 
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that of the Earth. Nevertheless, it is larger than M Pjm i g for this low- A;, case. Most cores migrate 
toward their host stars with M p less than a few Mq. Although these masses are < -Mp )trunc , gas 
accretion is too slow (eq. [11]) for them to acquire any significant amount of mass prior to gas 
depletion. The detection of close-in Neptune-mass around M dwarf stars and the modest detection 
frequency of gas giants around solar-type stars are inconsistent with the results of series B in Figure 
9a and therefore, we suggest that A v is substantially larger than unity. 

An upper limit on the magnitude of A v may be inferred from models xC. With A v = 100, 
the condition for the onset of migration becomes much more stringent and most cores do not have 
sufficient mass to undergo migration. For solar-type stars, a few relatively massive cores can form 
rapidly in disks with very large /d- These systems can migrate to form close-in Jupiter-mass gas 
giants. The mass distribution of close-in planets around solar type stars (model 1.0C) is skewed 
to 10 3 M e with lower cut-off below ~ 0.5Mj, which is inconsistent with the observed mass-period 
distribution of extra solar planets (Figure 1). For low-mass host stars (models 0.2C and 0.4C), 
many intermediate-mass planets can form during and after gas depletion. But they retain their 
initial semi major axis. The mass distribution of the close-in planets peaks near the mass of Saturn 
and hardly any planets have masses comparable to that of Neptune. These simulation results are 
again inconsistent with the observed mass-period distribution of extra solar planets. Furthermore, 
in series C, the deficit of planets with intermediate masses and periods is too pronounced to be 
consistent with observed one. Therefore, we infer A v ~ 10. 



6. Dependence on the disk mass 

The expressions in equations (1) and (2) indicate that the growth rate and asymptotic mass 
of cores are increasing function of Ed- In the standard series of models, we set h& = (M*/M & ) 2 . 
With this prescription, Ed of disks around low-mass stars is relatively small. Consequently, the 
emergence of gas giants occurs preferentially around massive stars. 

The dependence of the disks' Ed on the M* of their host stars is poorly known. On the 
theoretical side, gravitational instability may limit the amount of mass which can be retained by 
the disks, especially those around low-mass stars. But, smaller L^'s and lower intensity of ionizing 
photons may also reduce the influence of the magneto-rotational instability (Gammie 1996) and the 
angular momentum transfer efficiency so that more mass may be stored in disks around low-mass 
stars. Best available observational data suggest M oc M% (see §2.2) but the dependence of Ed on 
M* is poorly known. In view of these uncertainties, we introduce another two series of models. 

The parameters of models xD and xE are identical to those of models x, respectively. The 
only difference is that we set /id = M*/M Q in series D and it is set to be unity for all M* in series 
E. In comparison with the standard models, disks around low-mass stars are less deficient in these 
new models while those around the solar type stars remain the same. 

The predicted distributions for series D and E are shown in Figure 10. Since M Pim i g and M P)trunc 
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do not depend on nor S g , planets' migration and gas accretion undergo the same paths as in 
standard models for the same mass planets. As a result, the mass distribution of close-in planets is 
self-similar among models x, xD, and xE. Since massive disks exist around lower mass stars more 
frequently in series D (and even more frequently in series E), the amplitude of mass distributions 
around lower mass stars are enhanced in these models. In series D and E, inferred frequency of 
Jupiter mass planets around M stars are comparable to those around G stars (Figure 10). The 
sparse detection of close-in Jupiter-mass and Neptune-mass planets around M stars suggests that 
the disk masses are rapidly increasing function of M* as we have assumed in the standard models 
with h d = (M*/M ) 2 . 

7. Summary and discussions 

Observational discovery of extrasolar planets is advancing rapidly. We now have sufficient 
amount of data to carry out statistical characterization of planetary properties and to place con- 
straints not only on the dominant mode of planet formation but also the range of physical quantities 
which determine their growth and migration rates. In this paper, we focused our discussion on the 
mass function of close-in planets around stars with various masses, in particular lower masses than 
the solar mass, because its origin is determined by the delicate balance of various processes and 
they are the most conspicuous companions of nearby stars. Although many (> 90%) planets which 
once migrated to the proximity of their host stars may be eliminated (Paper II), we can compare 
the mass distribution of close-in planets among around stars with various masses. If the elimination 
factor is taken into account, rough comparison is also possible between close-in and more distant 
planets. 

Our results are summarized as follows. 

1. Dynamically-isolated, close-in, Neptune- mass planets with silicate and iron cores can form 
in relatively massive disks around solar-type stars. But their frequency is expected to be an 
order of magnitude smaller than that of close-in Jupiter-mass planets. 

2. Dynamically-isolated, close-in, Neptune-mass ice giants can form in lower-mass stars. Their 
frequency peaks around the M dwarfs. Since the luminosity of M dwarfs is weak, the ice 
boundary is located well inside 1 AU. These planets are formed at around 1AU but outside 
the ice boundary and their cores are primarily composed of volatile ices. Planets which 
migrated to the stellar proximity with masses in the range of 5 — 15M® may acquire, in 
situ, a limited amount of additional gas, but gas accretion is immediately quenched by gap 
formation because of small aspect ratio of the disk in the proximity of the host star. Because 
these planets compose mostly of icy material and M stars' luminosity is relatively weak, they 
may have water- vapor atmosphere and water ocean (Leger et al. 2004). 

3. Embryos with mass lower than 10M e cannot migrate to the proximity of their F, G, and 
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K dwarf host stars through type II migration. Detection of dynamically-isolated Earth- 
mass close-in planets may be attributed to type I migration of low-mass embryos (Ward 
1997a), sweeping secular resonances (Lin et al. 2004), or sweeping mean motion resonances. 
Around late M dwarfs, however, dynamically-isolated, a few earth-mass planets can form with 
temperature comparable to that of the earth. 

4. Around M dwarfs, the formation probability of gas giants is much reduced. The relatively 
low Ed prevents the emergence of sufficiently massive cores prior to the severe depletion of 
gas in their nascent disks, which is consistent with the results by Laughlin, Bodenheimer & 
Adams (2004). 

5. The mass function of close-in planets generally have two peaks at about Neptune mass and at 
about Jupiter mass. The lower-mass peak takes the maximum frequency for M stars, while 
the higher-mass peak is much more pronounced around higher-mass stars (F, G, K dwarfs). 
These are because planets tend to undergo type II migration after fully accreting gas around 
the higher-mass stars while they tend to migrate faster than gas accretion around M stars. 
Unless the termination location of planetary migration is a decreasing function of M*, close-in 
Neptune-mass planets around M dwarfs are easier to detect than those around G dwarfs. 

6. The mass function of dynamically- isolated close-in planets around stars with various masses 
can also be used to calibrate the sufficient condition for the onset of planetary migration and 
for the termination of gas accretion due to planet-disk tidal interaction. 

The metallicity dependence on frequency of extrasolar gas giant planets may not be easily 
accounted for by the gravitational instability scenario, e.g., (Boss 2001), while it is naturally ac- 
counted for by the sequential core accretion scenario that we are based on (Paper II). The condition 
of the gravitational instability is 1 > Q = CsQ.-^/ttGY^^ oc m}^ 2 /E g = M* - ^, and its radial 
wavelength A = 27r 2 G£g/f^ oc M*~ l (Toomre 1964) where E g oc M* . Hence, the instability may 
be more limited and result in smaller clumps around lower- mass stars, which could account for the 
above features 4 and 5, if (3 = 2 is assumed as in the standard series in the present paper. However, 
if (3 = 1, the gravitational instability scenario cannot account for the above features 4 and 5, while 
the sequential core accretion scenario still shows the tendency for the features. More detailed study 
on the dependence of disk mass on stellar mass M* is needed. 

Our work is primarily motivated by the discovery of GJ436b (Butler et al. 2004). Our theo- 
retical extrapolations can be tested with the following statistical properties of close-in planets to 
be discovered by various techniques. 

1. As mentioned above, the mass function of dynamically- isolated close-in planets as a function 
of the spectral classes of their host stars is particularly useful in the determination of the 
growth, migration, and disk depletion time scales. 
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2. A comparison between the frequencies of gas giants with close-in orbits and those with ex- 
tended orbits can provide constraints on the migration condition, survival criteria, and disk 
mass as functions of the host stars' mass. 

3. In systems with multiple giant planets, the dynamical architecture may provide clues on 
whether the migration of the close-in planets are driven by planet-disk tidal interaction or 
sweeping secular resonance. 

4. A comparison of atmospheric properties of close-in Neptune-mass planets around G dwarfs to 
that around M dwarfs can verify the conjecture that the former have silicate and iron cores 
whereas the latter have ice cores. 

The following observations of protostellar disks may provide useful input to the model: 

1. A spatially resolved image of disks can directly provide information of £d and the temperature 
distribution around any given host star. 

2. The dependence of S<j on the mass of the host stars determines the functional form of hd- 

3. A relation between the disk mass and accretion rate onto the host stars places a constraint 
on the rate of type II migration. 

4. A direct measurement of the gas distribution is particularly important in determining Td ep 
and h g . 

On the modeling side, we need to consider: 

1. The possibility of radiative feedback on the termination of gas accretion. 

2. The stoppage of type II migration and the survival of short-period planets. 

3. The rate (and direction) of type I migration of cores. 

4. The enhanced probability of multiple planet formation. 

5. Effect of dynamical interaction between multiple planets during and after gas depletion. 

6. The radial distributions of and S g . Effects of different power-law index of a dependence 
from —1.5 that we used. More realistic time evolution of S g . 

7. The influence of a stellar companion on the emergence and survival of planets. 
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Some of these issues will be addressed in future discussions. 

We thank G. Laughlin, M. Nagasawa, G. Ogilvie, and S. Vogt for useful discussions, and the 
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under its Origins program, JPL 1228184 under its SIM program, and NSF through AST-9987417 
and by JSPS. 



-24- 



REFERENCES 

Aarseth, S. J., Lin, D. N. C, Palmer, P. L. 1993, ApJ, 403,351 
Beckwith, S. V. W., & Sargent, A. L, 1996, Nature, 383, 139 
Balmforth, N. J. k Korycansky, D. G. 2001, MNRAS, 326, 833 
Bodenheimer, P., & Pollack, J. B. 1986, Icarus, 67, 391 
Bodenheimer, P., Hubickyj, O. & Lissauer, J. J. 2000, Icarus, 143, 2 
Boss, A. P. 2001, ApJ, 551, L167 

Boss, A. P., Wetherill, G. W. & Haghighipour, N., 2002, Icarus, 156, 291 

Bryden, G., Xingming, C, Lin, D. N. C, Nelson, R., & Papaloizou, J. C. B. 1999, ApJ, 514, 344 
Butler, P. et al. 2004, ApJ, in press 

D'Alessio, P, Calvet, N., Hartmann, L. 2001, ApJ, 553, 321 
Goldreich, P., & Tremaine, S. 1980, ApJ, 241, 425 

Greenberg, R., Hartmann, W. K., Chapman, C. R., & Wacker, J. F. 1978, Icarus 35, 1 
Gammie, C. F. 1996, ApJ, 457, 355 

Gu, P., Lin, D. N. C, & Bodenheimer, P. H. 2003, ApJ, 588, 509 
Haisch, K. E., Lada, E. A. & Lada, C. J. 2001, ApJ, 553, L153 
Hayashi, C. 1981, Prog. Theor. Phys. Suppl., 70, 35 

Hayashi, C, Nakazawa, K. k Nakagawa, Y. 1985, Protostars and Planets II, ed. D. C. Black k 
M. S.Matthew (Tuscon: Univ. of Arizona Press), 1100 

Ida, S., k Makino, J. 1993, Icarus, 106, 210 

Ida, S., Bryden, G., Lin, D. N. C. k Tanaka, H. 2000, ApJ, 534, 428 

Ida, S. k Lin, D. N. C. 2004, ApJ, 604, 388 (Paper I) 

Ida, S. k Lin, D. N. C. 2004, ApJ, 616, 567 (Paper II) 

Ikoma, M., Nakazawa, K. k Emori, E. 2000, ApJ, 537, 1013 

Ivanov, P. B., Papaloizou, J. C. B. k Polnarev, A. G. 1999, MNRAS, 307, 791 

Iwasaki, K., Emori, H., Nakazawa, K. k Tanaka, H. 2002, PAS J, 54, 471 



-25- 



Kokubo, E. k Ida, S. 1996, Icarus, 123, 180 

— . 1998, Icarus, 131, 171 

-. 2000, Icarus, 143, 15 

-. 2002, ApJ, 581, 666 
Roller, J., Li, H. k Lin, D. N. C. 2003, ApJ, 596, L91 
Kominami, J. k Ida, S. 2002, Icarus, 157, 43 
Laughlin, G., Steinacker, A. k Adams, F. C. 2004, ApJ, 608, 489 
Laughlin, G., Bodenheimer, P. k Adams, F. C. 2004, ApJ, 612, L73 
Leger et al. 2004, Icarus, 169, 499 

Lin, D. N. C. k Papaloizou, J. C. B. 1980, MNRAS, 191, 37 

Lin, D. N. C. k Papaloizou, J. C. B. 1985, Protostars and Planets II, ed. D. C. Black k 
M. S.Matthew (Tuscon: Univ. of Arizona Press), 981 

Lin, D. N. C, k Papaloizou, J. C. B. 1993, in Protostars and Planets III, ed. E. H. Levy and J. I. 
Lunine (TucsomUniv. of Arizona Press), 749 

Lin, D. N. C, Bodenheimer, P. k Richardson, D. 1996, Nature, 380, 606 

Lin, D. N. C, Nagasawa, M. k Thommes, E. W., submitted 

Lissauer, J. J. 1987, Icarus, 69, 249 

Lissauer, J. J. k Stewart, G. R. 1993 

Lubow, S. H., Seibert, M. k Artymowicz, P. 1999, ApJ, 526, 1001 

Marcy, G. W., Cochran, W .D., k Mayor, M. in Protostars and Planets IV, ed. V. Mannings, A. P. 
Boss and S. S. Russell (TucsomUniv. of Arizona Press), 1285 

Malhotra, R. 1993, Nature, 365, 819 

Mayor, M. k Queloz, D. 1995, Nature, 378, 355 

McArthur, B. 2004, ApJL, 614, L81 

Mizuno, H. 1980, Prog. Theor. Phys. Suppl., 64, 54 

Mohanty, S., Jayawardhana, R. k Navascues, 2003, ApJ, 593, L109 

Muzerolle, J. et al., 2003, ApJ, 592, 266 



-26- 



Nagasawa, M., Tanaka, H. & Ida, S. 2000, AJ, 119. 1480 
Natta, A. et al., 2004, A&A, 424, 603 

Nelson, R. P., Papaloizou, J. C. B., Masset, F. & Kley, W. 2000, MNRAS, 318, 18 
Nelson, R. P. & Papaloizou, J. C. B. 2004, MNRAS, 350, 849 
Papaloizou, J. C. B. & Lin, D. N. C. 1995, ARA&A, 33, 505 

Pollack, J. B., Hollenbach, D., Beckwith, S., Simonelli, D. P., Roush, T., & Fong, W. 1994, ApJ, 
421, 615 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & Greenzweig, Y. 1996, 
Icarus, 124, 62 

Rafikov, R. R. 2003, AJ, 125, 922 

Sano, T. et al. 2000, ApJ, 543, 486 

Santos, N. et al. 2004, A&A, 426, L19 

Safronov, V. 1969, Evolution of the Protoplanetary Cloud and Formation of the Earth and Planets 
(Moscow: Nauka Press) 

Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 

Stevenson, D. J. 1982, P&SS, 30, 755 

Tanaka, H. & Ida, S. 1997, Icarus, 125, 302 

Tanaka, H., Himeno, Y., & Ida, S. 2004, submitted 

Thi, W. F., Blake, G. A., van Dishoeck, E. F., van Zadelhoff, G. J., Horn, J. M. M., Becklin, E. 
E., Mannings, V., Sargent, A. I., van den Ancker, M. E.; Natta, A. 2001, Nature, 409, 60 

Toomre, A. ApJ, 139, 1217 

Ward, W. 1981, Icarus, 47, 234 

Ward, W. 1986, Icarus, 67, 164 

Ward, W. 1997, Icarus, 126, 261 

Ward, W. 1997, ApJ, 482, L211 

Wetherill, G. W. 1980, ARAA, 77, 330 

Wetherill, G. W. & Stewart, G. R. 1989, Icarus, 77, 330 



-27- 



Wuchterl, G., Guillot, T. & Lissauer, J. J. 2000, Protostars and Planets II, ed. V. Mannings, A. P. 
Boss k, S. S. Russell (Tuscon: Univ. of Arizona Press), 1081 

Wyatt, M. C. et al. 2003, MNRAS, 342, 876 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



-28- 



10' 
10 E 



1 1 1 II Ml 

! all " [ 


1 1 J 

° • ■ v.° . 

























10" 



10" 



10 1 



m 

CO 



cd 

CP 



i i rT 

! M 



I 0" | 

I 10 2 | 

i in' i 



10 J 

io 2 k 
io 1 r 
l - 



io 3 k 

10 2 f 

io 1 t 



10 ; 
10 J 
IO 1 



K 



! G 



10 



10" 



IO 1 



a [AU] 



Fig. 1. — The distributions of semimajor axis (a) and mass (M p sku) of discovered extrasolar 
planets. Unit of mass is Earth mass M e (Jupiter mass is Mj ~ 320M e ). The data are taken from 
"The Extrasolar Planets Encyclopedia" (http://cfa-www.harvard.edu/planets/) as of February, 
2005. The planets around subgiants and those discovered by transit survey are excluded. The right 
panels show the planets around M, K, G, and F stars, respectively. The left panel shows all the 
data. The dotted lines show observational limits (v r = lOm/s) for doppler survey. Detection of 
larger a (> 3AU) planets, which have longer orbital periods, is also limited. 
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Fig. 2. — The distributions of semimajor axis (a) and mass (M p ) of planets predicted by the Monte 
Carlo simulations, (a) final semimajor axes (ag n ) and masses (M p ^ n ) at t = 10 9 yrs, and (b) initial 
semimajor axes (ami) and M Pjmi . The symbols x (= 0.2,0.4,0.6, 1.0 and 1.5) represent the host star 
mass scaled by the solar mass, M*/Mq. In (b), close-in planets with a^ n < 0.05AU are marked by 
black crosses, while the other planets are marked by gray dots. The thin solid black lines indicate 
the isolation mass M Ci ; so with fa = 30. The thick solid black lines express the critical mass for 
radial migration, M Pirmg = A u M^y ls with A u = 10. The dashed lines are the truncation mass for 
gas accretion, M Pj t rU nc = ^4th-^g,th with A t h = 3.4. 
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(a) a fm <0.05AU (b) 0. 1 AU<a fm < 1 AU 




M pfin [earth mass] M p (in [earth mass] 



Fig. 3. — The calculated distribution of final mass of planets for (a) close-in planets at ag n < 0.05AU 
and (b) planets at 0.1AU < aa n < 1AU (filled circles). The labels x are the same as Figure 2. Since 
it is expected that most of close-in planets may fall onto their host stars, the calculated amplitude 
iV in (a) is reduced by a factor 10. Observed distributions are also plotted with open triangles. The 
sin i factor is neglected for simplicity (it enhances the observed values only by A/n on average). The 
number of runs in each model does not reflect the number of targets for current doppler survey. 
Also note that observed distributions do not exactly correspond to host stars' mass of each model. 
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Fig. 4. — The fdhd distribution we used for the standard model, which is a gaussian distribution 
in terms of log 10 fd with a center at log 10 fd = and dispersion of 1. In the standard model, we 
assume h d oc (M*/M ) 2 . We omit the high h d f d tail at > 3O(M*/M ), since such heavy disks 
are self gravitationally unstable. Filled circles, open circles, filled squares, open squares and filled 
triangles represent the cases of M* = 0.2, 0.4, 0.6, 1.0, and 1.5M . 
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Fig. 5. — The mean mass and the characteristic mass associated with the peak of the mass dis- 
tribution given in Figures 3, as a function of host stars' mass M*. (a) the close-in planets with 
a fin < 0.05AU and (b) the planets at 0.1 AU < ag n < 1AU. The mean mass and the peak mass 
are plotted with filled circles and open squares. For (b), only planets with M p over a deficit 
(M p > 50M®) are considered (see Figures 2a and 3b). 
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Fig. 6. — An example of formation of a close-in Neptune-mass planet around an M star. A plant is 
initially at 0.74AU in a disk with faka = 7.2 and Td cp = 9.2 Myr around a star with M* = 0.4M Q . 
The upper and middle panels show time evolution of mass and semimajor axis of the planet. 
In the upper panel, the isolation mass of a core (M Ci ; so ), its critical mass for initiation of gas 
accretion (M Ci h y dro)> the planet's critical mass for migration (M g;m i g ) and the truncation mass of 
gas accretion (M git runc) are plotted for comparison. In the lower panel, core accretion time scale 
(7"c,grow = M c /M c ), gas accretion time scale (tkh) and migration time scale (r m i g ) are plotted. 
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Fig. 7.- 
a K star, 
with M* 



An example of formation of an intermediate-mass planet with intermediate period around 
A plant is initially at 1.8AU in a disk with f^h^ = 2.0 and Td ep = 2.3 Myr around a star 
= 0.6M©. Notations are the same as Figure 5. 
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Fig. 8. — The theoretically predicted fraction (r/j) of stars which bears giant planets currently 
detectable with current doppler survey, as a function of their mass (M*). The results of the 
standard model with = (M*/M & ) 2 are plotted by filled circles, while the results of series D and 
E with /id = M^/M & and 1 are plotted by open squares and triangles, respectively. The detectable 
conditions are radial velocity v r > lOm/s and period P < 4 years. We do not include planets 
stalled at 0.04AU in the evaluation of r/j, since most of them may fall onto their host stars. The 
distribution of /a is the same as in Figure 4. For each /d, a is selected as \og(a,j + i/aj) = 0.2 

(i = 1,2,...). 
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Fig. 9. — The predicted distributions of semimajor axis (a) and mass (M p ) of planets in the final 
state at t = 10 9 yrs: (a) model xB with M Pim i g = M gjV i s , (b) model xC with M Pim i g = 100M giV i s , 
where x in xB and xC represent the results with M* = xM Q . 
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Fig. 10. — The same as Figure 9 except (a) model xD with = M*/M & , (b) model xE with 
h d = l. 



